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We report an ultrasonic investigation of the elastic moduli on a single crystal of hexagonal YMnOa 
as a function of temperature. Stiffening anomalies in the antiferromagnetic Neel state below Tjv = 
72.4 K are observed on all the four elastic moduli Cu. The anomalies are the most important on 
Cii and Cee for in-plane elastic deformations; this is consistent with a strong coupling of the lattice 
with the in-plane exchange interactions. We use a Landau free energy model to account for these 
elastic anomalies. We derive an expression which relates the temperature profile of the anomaly to 
the order parameter; the critical exponent associated to this parameter /3 — 0.42 is not consistent 
with a chiral XY or 3D Heisenberg universality class, but more in agreement with a conventional 
antiferromagnetic long range order. A tiny softening anomaly on Cu for which hysteresis effects are 
observed could be indicative of an interaction between ferroelectric and magnetic domains at Tat. 
Moreover, magnetic fluctuations effects both above and below Tm are identified through abnormal 
temperature and magnetic field effects. 

PACS numbers: 75.47.Np, 75.50.Ee, 75.40.Cx 



I. INTRODUCTION 

Hexagonal YMnOa is one of the most well-known ex- 
amples of multiferroics, a material that possesses both 
ferroelectricity and magnetism. The ferroelectric (FE) 
transition occurs at a relatively high temperature above 
900 K. Below the Curie temperature the compound crys- 
tallizes in a hexagonal structure with the space group 
P63CTO The Mn'^'*' ions form nearly triangu- 

lar networks in z—0 and z—1/2 layers which stack in the 
ABAB sequence along the z axis with a wide separation 
introduced by the intervening Y and O ions This 
suggests a predominant two-dimensional (2D) character 
in the ab plane. For this system the geometrical frustra- 
tion of the antiferromagnetic (AF) spins on the triangular 
lattice results in the 120° spin ordering below the Neel 
temperature Tn '--^ 70 K. 

The Curie-Weiss (CW) temperature deduced from 
susceptibility measurements is approximately ten times 
higher than Tn 0, S H and the magnetic moment is 
lower than the expected value from the fully polarized 
This reduction in Tjv and magnetic mo- 
ment is consistent with strong spin fluctuations due to 
geometrical frustration and/or low dimensionality. The 
coexistence of novel magnetism and ferroelectricity sug- 
gests a peculiar interaction between the spins and the lat- 
tice degrees of freedom, although the low Tjm, compared 
to the ferroelectric transition, infers weak coupling be- 
tween magnetism and ferroelectricity. Although various 
crystallographic, optical and magnetic structural studies 
have been performed to characterize independently the 
FE and magnetic properties P, d, @], only a few ones 
have been directed to understand their interrelation. Di- 
electric constant Q , specific heat 0, 0] and thermal con- 



ductivity 0, [T3] experiments revealed a spin-lattice cou- 
pling whose understanding remains an open question. As 
far as magnetism is concerned, neutron scattering studies 
[ill fl^ Il3l . Il4 | have reported the presence of unconven- 
tional spin fluctuations both above and below Tn that 
could be related to the anomalous increase of the thermal 
conductivity upon magnetic ordering [9^]. However, the 
question of coexistence of long-range order with signifi- 
cant spin-fluctuations in the Neel state for the triangular 
lattice remains. The symmetry of the order parameter 
and the nature of the universality class associated with 
the 120° spin structure are also not resolved up to now. 

To investigate the coupling between magnetism and 
the lattice degrees of freedom in YMnOa, we report an 
ultrasonic investigation of the elastic moduli of a single 
crystal in the temperature range 2-150 K, both above and 
below the spin ordering temperature Tn- Although elas- 
tic anomalies at Tn are observed on all studied moduli, 
an important coupling of the magnetic order parameter 
is only obtained with an in-plane strain field. A Landau 
free energy expansion is used to interpret these anoma- 
lies and a critical analysis of the temperature behavior 
just below Tn is tentatively given. Small temperature 
hysteretic effects are observed indicating a domain struc- 
ture. A peculiar temperature profile above Tn for two 
moduli is discussed in relation with the magnetic fluctu- 
ation issue. 



II. EXPERIMENT 

High quality YMnOs single crystals are grown using 
the standard floating zone technique. Powder X-ray 
diffraction measurements confirmed that the crystal in- 
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vestigated here has the PGscm hexagonal structure at 
room temperature (a=&=6.1380(l)A, c=11.4045(l)A). A 
single crystal was then oriented with a back Laue X- 
ray diffraction technique and cut with a wire-saw to ob- 
tain two sets of parallel faces respectively oriented per- 
pendicular (separation distance 3.29 mm) and parallel 
(3.54 mm) to the hexagonal c axis. The faces were sub- 
sequently polished to obtain a mirror like aspect. We 
used a pulsed-echo ultrasonic technique in the reflection 
mode to measure the velocity and amplitude of longi- 
tudinal and transverses waves propagating along both 
directions. The plane acoustic waves were generated 
with LiNbOa piezoelectric transducers (fundamental fre- 
quency 30 MHz and odd overtones) bonded to one paral- 
lel face with silicon seal. The technique consists in mea- 
suring the phase shift and the amplitude of the first elas- 
tic pulse reflected from the opposite parallel face. Doing 
so, this technique yields only velocity and attenuation 
variations; an absolute value of the velocity is obtained 
by measuring the transit time between different reflected 
echoes. The crystal was mounted in a Variable Tem- 
perature Insert ( VTI) ; the useful temperature range was 
2-150 K, the maximum value being imposed by phase 
transformation of the silicon seal around 200 K. A mag- 
netic field up to 14 Tesla could be applied along different 
directions. 

For a crystal having hexagonal symmetry the elastic 
modulus matrix has only five non-zero independent com- 
ponents Cij [l^. Four of these components can be ob- 
tained by propagating ultrasonic waves parallel and per- 
pendicular to the hexagonal c axis and by measuring the 
velocity of these waves with various polarizations (longi- 
tudinal and transverse) . For propagation along c, the lon- 
gitudinal and transverse velocities are respectively given 
by Vli — (Caa/p)^ and Vri = (6*44/^)2; for the perpen- 
dicular direction, besides the longitudinal one, two trans- 
verse velocities (polarization _L c and || c) are obtained: 
Vl2 = (Cii/p)^, Vt± = (Cee/p)' and Vt\\ = [C^i/p)^. 
By using the density p = 5.1 g/cra^ and the experimental 
velocities, the four Cij are obtained. We notice that Vti 
and Vt\\ are equal if the crystal is correctly oriented, a 
situation that was checked in our experiment. The abso- 
lute values of the moduli at low temperatures are given 
in Table I. 



TABLE I: Elastic moduli of YMnOs (units at 4K. 



Cii 


C33 


C44 




18.5 ±0.1 


29.8 ±0.5 


9.86 ±0.15 


5.94 ±0.05 




III. 


RESULTS 





The coupling between the lattice and the spin degrees 
of freedom generally yields anomalies on selected elastic 



moduli when three-dimensional magnetic order sets in. 
The selection is imposed by the nature of the coupling 
between the elastic deformation and the order parameter, 
coupling which is highly sensitive on the symmetry of the 
latter. This coupling is confirmed for all elastic moduli 
investigated in the present work on multiferroic hexago- 
nal YMnOa. Our ultrasonic velocity measurements have 
been performed at different frequencies between 30 and 
300 MHz; no significative frequency dependence were no- 
ticed on the velocity. The data presented below were all 
obtained at 103 MHz. 
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FIG. 1: Temperature profile of AC33/C33 and AC44/C44. Tn 
is indicated by an arrow. 

In figure 1 we present the temperature profile of the rel- 
ative variation of the two elastic moduli that are the least 
affected by magnetic ordering, ACaa/Caa and AC44/C44. 
When the temperature is decreased from 150 K, the mod- 
ulus Caa which implies a pure compression along the c 
axis smoothly increases and saturates at low tempera- 
ture. This temperature profile appears normal in view 
of the expected stiffening following the progressive dis- 
appearance of phonon anharmonic effects at low temper- 
atures. There is however a small increase of slope below 
72 K indicating a further stiffening when AF magnetic 
ordering occurs at T/v- CI44, related to a shear deforma- 
tion in a plane containing the c axis, presents a similar 
temperature profile. There is however a much more pro- 
nounced upturn of the modulus below 72.5 K confirming 
a definite stiffening below Tm 0.1%). 
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FIG. 2: Temperature profile of ACn/Cn and ACee/Ce 



The temperature profile of the relative variation of the 
other two moduli appears both similar and different from 
the previous ones and, in a sense, less normal. These 
data are presented in figure 2 over the same tempera- 
ture range. Cn, corresponding to a pure compression 
in the ab plane, increases weakly below 150 K and, on 
this scale, appears to saturate around 90 K. Then, an 
abrupt stiffening begins at T/v = 72.4 K until saturation 
is obtained below 15 K; the amplitude of the stiffening 
is approximatively one percent. Cgg is also associated to 
in-plane deformations and it shows a non conventional 
behavior. The modulus decreases below 150 K until a 
minimum is reached around 75 K; then it increases fur- 
ther (~ 3%) with a clear acceleration near 72 K. Below 
72 K the thermal profile resembles the one observed on 
Cii. 

It appears then that the four elastic moduli investi- 
gated here, Cn, C33, C44 and Cge stiffen progressively 
below the Neel temperature Tn = 12A K. Moreover, 
the temperature profile of the anomalies seems to mimic 
the AF rnagnetic order parameter, the magnetic moment 
m(r) I, [ll| . It is worth to notice that significant elastic 
anomalies are observed for moduli implying elastic de- 
formation (compression and shear) only in the ah plane, 
Cn and Cge; these are also the only moduli to present an 
unconventional elastic behavior above T/y. These obser- 
vations could appear consistent with a two-dimensional 
character of the magnetic properties. Moreover, since the 



anomaly on C33 is very small, this means that inter-plane 
exchange interactions and are not strongly coupled 
to the lattice compared to in-plane ones Ji and J2 which 
yield an important anomaly on Cn in agreement with 
the analysis of spin fiuctuations 12]. The analysis of 
the temperature dependence of these anomalies is highly 
delicate since it necessitates a precise knowledge of the 
overall temperature background both above and below 
T/v. Such an analysis appears difficult for C33 because 
the amplitude of the anomaly is too small compared to 
the total modulus variation with temperature. However 
this could be done for Cn and Cee since the amplitude 
of the anomalies, respectively around 1% and 3%, are 
much larger. It will also be tempted for C44 even if the 
variation is only around 0.1%. 
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FIG. 3: Temperature profile of the relative variation of the 
elastic modulus Cn: (a) reduced scales for T > Tat; (b) hys- 
teresis near Tm, arrows indicating the direction of the tem- 
perature sweeps. 

Before doing so, let us examine more precisely the tem- 
perature profile of Cn just above T/v. The relative vari- 
ation ACn/Cn (relative to the value at 100 K) is pre- 
sented in figure 3a. The modulus Cn increases smoothly 
as the temperature is decreased from 100 K; then it goes 
through a maximum around 82 K, reaches a minimum 
just above 75 K and increases further before reaching the 
AF transition. A sharp downward step (0.02%) occurs 
below 72.8 K before an abrupt upturn at 72.4 K signaling 
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a progressive stiffening. These features in this tempera- 
ture range were observed at all frequencies without any 
significant modification. This abrupt downward step pre- 
ceding a stiffening is accompanied by a small hysteresis 
(0.4 K) when sweeping the temperature up and down 
through the transition as shown in figure 3b. The tem- 
perature profile above T^r is only modified up to 75 K 
where a local minimum is observed (Fig. 3a). Although 
they show also hysteresis, the other moduli do not show 
a downward step at Tjq . A minimum is however observed 
on Cgg near 75 K (Fig. 2), minimmn that could be a sig- 
nature of Cii since Cge — \{Cii — C12). 
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FIG. 4: Relative variation of the elastic modulus Cu, C44 
and Cee as a function of the reduced temperature (1 — T/Tjv). 
Straight lines represent fits to Eq.(7). 



IV. DISCUSSION 

The effects of magnetic frustration in the stacked tri- 
angular lattice have been the object of many studies over 
the years [3 and ultrasonic techniques [13, [H, [3| have 
been used extensively as a mean to obtain the magnetic 
phase diagram of ABX3 systems. Similarly to what has 
been done in these systems, we can use a phenomeno- 
logical Landau model to take account of the observed 
changes in the elastic properties of YMnOs. If one takes 
into account the coupling between the order parameter 
Q and the strain e^, three terms contribute to the free 
energy F(Q,ej). 

F{Q,e,) = L{Q) + F,i{e,) + F,{Q,e,) (1) 

where L{Q) is the Landau-type free energy expressed in 
terms of a power series of the order parameter Q, Fei{ei) 
the elastic energy associated to the deformations and 
Fc{Q, Si) the coupling between the strain and the order 
parameter Q. Because of the 120° spin arrangement in 
the ab plane [1, [13 in the absence of a magnetic field, we 



can write with in-plane spin component S±^ as the order 
parameter (we neglect any Sz component) , 

L{SA_) = aSl+hS\ (2) 

If we neglect any terms associated to the elastic moduli 
C12 and Ci3jthe elastic energy for a hexagonal structure 
is given by 

Fei{e,) =i Cii(e? + el) + \ C^^el + \ C^^{el e^) 

+ 5^666^ (3) 

where Cy represent the bare elastic moduli in the high 
temperature phase. For the coupling term, time inversion 
symmetry imposes a quadratic dependence on the order 
parameter, S\: the coupling term associated to the strain 
Ci can thus be written as 

F,{S^,e,)^g,^,Sle\ (4) 

where gi^r is the coupling constant associated to the de- 
formation Ci elevated to the power r. The term with 
r = 2 is always allowed whatever the symmetry but r = 1 
implies inequivalence of |ei| and — |ei|. If we minimize the 
Landau energy relative to the strain e^, we obtain a spon- 
taneous strain ei(S\}) which is used to derive the elastic 
constant ■ The biquadratic coupling term (r = 2) 

leads to a stiffening of the elastic constant Ca according 
to the general form 

C; = a. + 2g,.2S\ (5) 

For the case of linear-quadratic coupling (r = 1), a neg- 
ative step is obtained according to 

= C^^ - 1^ (6) 

The results presented in figures 1 and 2 are consistent 
with the prediction of Eq.(5), a stiffening related to . 
We give in Table II the relative value of the coupling 
constant 3^,2 as deduced from figures 1 and 2. The largest 
constant is obtained for C\\ and g\^i is thus taken as 1; 
56,2 associated to Cge, is also very important around 0.9. 
The constants 53,2 and 54^2 are found smaller by a factor 
between 20-25. This anisotropy of the coupling to the 
order parameter appears consistent with the one observed 
in the in-plane and inter-plane exchange constants (l^ 
and with a XY symmetry. 



TABLE II: Relative values of the coupling constant 3^,2 • 





53,2 


34,2 


96.2 


1.0 


0.04(1) 


0.05(1) 


0.90(5) 



We have then tried to analyze the variation of the elas- 
tic moduli Cii, C44 and Cge relative to their value at T/v 
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by fitting the order parameter to S± ^ {1 — T/T^Y near 
the critical temperature T/y. 

^ = D,{l-T/T^f^ (7) 

where Di is a constant depending on gi^ijCa. This analy- 
sis for Cii and Cge is surely more valid than for C44 since 
the anomalies are at least 10 times larger at T/v. The 
results of this analysis are presented in figure 4 on a log- 
log scale. Not only the overall temperature dependence 
appears very similar for Cn and Cgg, but the same expo- 
nent 2/3 = 0.85(5) is obtained with Tn = 72.43(5). Sur- 
prisingly the same exponent is also found for C44, what 
means that the overall temperature profile of the back- 
ground near T/v is not important to determine the critical 
behavior. Also hysteresis effects were found to only af- 
fect Tn (Fig.3), not the exponent. The value /3 = 0.42(3), 
which is obtained over only two decades of the reduced 
temperature, is neither consistent with chiral XY (0.25) 
and chiral Heisenberg (0.30) universality classes or with 
XY (0.35) and Heisenberg (0.36) ones Although this 
appears to contradict our previous observation pertain- 
ing to the coupling constants, it is clear from Fig. 4 that 
the same physical phenomenon is at the origin of the elas- 
tic anomalies on C\\ and Cge, and by extension on the 
other two moduli. The AF magnetic order parameter 
S\_ has been associated with 3D Heisenberg or 3D XY 
symmetry classes from specific heat j7| and muon-spin 
relaxation [2l| studies; our study is however not able to 
resolve this issue. 

Finally, the sharp downward step observed on C\\ 
(Fig. 3a) could be consistent with the prediction of Eq.(6). 
On the one hand, the observation of hysteresis in the 
vicinity of this peak could be consistent with weak fer- 
romagnetism, although this has been observed only in 
ScMnOs [3[ . On the other hand, coupHng between mag- 
netic and electric domains in YMnOa has been reported 
by imaging with optical second harmonic generation [2^ ; 
it has been suggested that the small anomaly observed 
on the in-plane dielectric constant [2^ could result from 
such coupling effects on domain walls [23|. The small 
peak on C\\ at T/v together with the observation of hys- 
teresis could result from an interaction between in-plane 
compression waves and the domain walls. 

Previous neutron scattering studies [H, [H, [13] re- 
ported the presence of unconventional spin fluctuations 
both above and below T/v. These fluctuations could ex- 
plain the anomalous behavior of C\\ (and by extension 
Cge) well above T^r as shown in figure 3a. Since similar ef- 
fects are not observed for C33 and C44, fluctuations pref- 
erentially 2D-like associated to in-plane exchange con- 
stants are likely responsible. The ultrasonic attenuation 
is generally more sensitive to fluctuations since it is a 
dissipative phenomenon. These attenuation variations 
Aa are difficult to obtain with precision on a wide tem- 
perature range and this is why they must be interpret 
with great care. We present in figure 5 the variation 
of the attenuation associated to the C\\ modulus as a 
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FIG. 5: Variation of the attenuation Aa of in-plane compres- 
sion waves a function of temperature at two frequencies, 103 
and 168 MHz. 

function of temperature at two frequencies. As the tem- 
perature is decreased from high temperature, the atten- 
uation generally decreases because the scattering of the 
acoustic waves by phonons is reduced. Here the atten- 
uation shows a minimum at a rather high temperature 
near 110 K for both frequencies and it increases for lower 
temperatures, a behavior that is not consistent with de- 
creasing phonon scattering. At 103 MHz, the AF mag- 
netic ordering produces a very sharp peak at Tjv but 
does not change significantly the monotonous increase of 
the attenuation that begun near 110 K and the attenu- 
ation increases down to the lowest temperature. At 168 
MHz, the temperature profile from 150 K to 72 K is con- 
served; however, the peak is more important at T/v and 
the attenuation decreases rapidly below the Neel tem- 
perature. At lower temperature, there is no saturation 
of the attenuation but a small increase below 30 K is 
observed. These attenuation data strongly suggest that 
spin degrees of freedom, in other words magnetic fluc- 
tuations both above (at least up to 110 K) and below 
Tjv, dissipate acoustic energy. At T/v, the sharp increase 
is due to the spontaneous strain yielding to eq.(6) and 
for which hysteresis is observed. Below Tjv, spin degrees 
of freedom are lost when 3D AF ordering occur and the 
attenuation decreases rapidly. The data reveal an in- 
crease of attenuation at temperatures much higher than 
T/v and a reduction below when 3D AF magnetic order- 
ing occurs. Although this reduction increases rapidly in 
amplitude with frequency, dissipation still occurs well be- 
low r^r. These observations are consistent with extraor- 
dinary spin-phonon interactions deduced from thermal 
conductivity data ^] . How these fluctuations are related 
to the domain structure remains an open question. The 
presence of domains and fluctuations near T/v could com- 
plicate the analysis of the critical behavior presented in 
Fig. 4 and affect the determination of the exponent (3. 
The magnetic character of the peculiar behavior of Cn 
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FIG. 6: Relative variation of Cn as a function of the magnetic 
field at two temperatures well below and above Tn- 

above Tjv can also be revealed by a magnetic field analy- 
sis. Up to 14 Tesla we did not find any measurable varia- 
tion of Tn whatever the orientation of the field relative to 
the c axis was. This probably implies that the amplitude 
of the magnetic order parameter below Tn is highly rigid 
relative to the magnetic field. Nevertheless some mag- 
netic field effects can be measured on Cn over the 2-150 
K temperature range. These are presented in figure 6 as 
Cii variations relative to the value at zero field. The data 
are presented as a function of at two temperatures, 10 
K well below T/v and 100 K well above for a field applied 
parallel to the c axis. The field profile is very similar for 
both temperatures and the absolute value has not even 
decreased by a factor 2 when the temperature has been 
multiplied by ten. These effects have thus to be related 
to the overall background: a magnetic field enhances the 
softening of Cn both above and below Tjv confirming the 
existence of important magnetic fluctuations in this geo- 
metrically frustrated magnet, fluctuations that are likely 



enhanced by a magnetic field. It has been suggested from 
heat capacity data [3] that a single Einstein contribution 
could adequately fit the low-temperature deviation, indi- 
cating that there is no evidence of a anomalous magnetic 
contribution for T < T/y. Our data do not support this 
suggestion. 



V. CONCLUSION 

In conclusion, our ultrasonic experiments have iden- 
tified an important coupling between the AF magnetic 
order parameter and the lattice in hexagonal YMnOa. 
This coupling yield important stiffening anomalies below 
Tjv only for in-plane elastic deformations, compression 
and shear. This implies that the coupling is established 
via the modulation of the in-plane exchange interactions. 
The critical exponent f3 for the AF order parameter near 
Tjv is found higher than the expected value from XY and 
3D Heisenberg universality classes. This situation favors 
a conventional antiferromagnetic long range ordering as 
suggested from other measurements. However, the data 
confirm important spin fluctuations over a wide temper- 
ature range above and below the Neel temperature. Fi- 
nally, a sharp step appearing on Cn at Tjv, for which 
hysteresis is observed, could be consistent with a weak 
coupling between ferroelectricity and magnetism at do- 
main walls. 
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